Instituto Superior Técnico 

Departamento de Matemática 

Secção de Álgebra e Análise 



ANALISE MATEMÁTICA IV 

FICHA 6 - SÉRIES DE FOURIER E 
MÉTODO DE SEPARAÇÃO DAS VARIÁVEIS 



(1) Determine o desenvolvimento em série de Fourier das seguintes funções: 

(a) f(x) = x para x e [—1,1]; 

(b) f(x) = x + 1 para x e [-1, 1]; 

(c) f(x) = cos 3 x para x E [— ir, ir]. 

Resolução: 

(a) A série de Fourier de f é da forma 

oo 

f( x ) = TT + y, ( a k cos(kTTx) + bk sm(kjTx)) 
fc=l 

Como f(x) é ímpar, todos os ak 's são 0. Quanto aos bk 's, são dados pela fórmula 

bk = xsm(kTTx)dx 



2 / xsm(kTTx)dx porque x sm(hivx) é par 
Jo 



I x cos(kirx) 



kn 
2(-l) fc 



f 1 COs(kTTx) 

+ / 7 -dx 

o Jo kir 



kn 
2(_i)fc+i 



+ 



kir 
Conclui-se que o desenvolvimento de Fourier de f para x e [— 1, 1] é 

°°^ 2(-l) fc+1 

f(x) = ) sm(kTrx) 

^— ' kn 
fc=i 

(b) A função x já foi desenvolvida em série pelo que nos resta desenvolver a função 
constante igual a 1 no intervalo [— 1, 1]. Por unicidade do desenvolvimento de Fourier, 
há uma única escolha possível para a^ s e bk 's tais que 

oo 

1 = — + 2, ( a k cos(kirx) + bk sm(kirx)) 
fe=i 

Claramente uma escolha possível é ao = 2 e a^ = bk = para k > 1. Por unicidade 
conclui-se então que o desenvolvimento de Fourier da função constante igual a 1 é 

2 

1 = - 

2 



e portanto o desenvolvimento de Fourier de f é 

2 J^ 2(— l) fc+1 
f( x ) = õ + Z_/ í sin(fcvra;) 



2 ^^ /CTT 

fe=l 
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(c) O desenvolvimento de Fourier de f é da forma 

OD 

ao 



f( x ) = 1T + 51 ( flfc cos ( kx ) + b k sin(kx)) 



fc=i 

Por unicidade do desenvolvimento de Fourier, qualquer desenvolvimento desta forma 
que se obtenha será o desenvolvimento de Fourier. Pela fórmula de DeMoivre tem-se: 

(cos x + i sin x) = cos(3a;) + i sin(3x) 
cos x — 3 cos x sin a; + í(3cos x sin x — sin x) = cos(3a;) + isin(3x) 

donde, igualando as partes reais, 

cos x — 3cosa;(l — cos x) = cos(3a;) 

4 cos x = 3 cos a; + cos(3x) 

s 3 1 , . 

cos x = - cosx + - cos(3a;) 
4 4 

Uma vez que este é um desenvolvimento da forma pretendida, conclui-se que é este 
o desenvolvimento de Fourier. Isto é tem-se bk = para todo o k > 1, ãk = para 

k ^ 1,3 eai = |, a 3 = \. 

D 

Comentário: Nas alíneas (b) e (c) do exercício anterior poder-se-ia também ter utilizado 
as fórmulas integrais para calcular os coeficientes a^ e bk mas esse processo seria muito 
mais trabalhoso, principalmente na alínea (c). <) 



(a) f{x) = 



(2) Determine o desenvolvimento em série de senos das seguintes funções: 

x se < x < 1 
1 se 1 < x < 2 ' 
(b) f(x) = cosx para x e [0, 2ir]. 

Resolução: 

(a) O desenvolvimento de f em série de senos no intervalo [0, 2] é da forma 



OD i 

f( x ) = ^òfcsin í 
fe=i ^ 



knx\ 
"2"/' 



Os coeficientes bk são dados pela fórmula 



bk = 



f(x) sin 



kirx 



(Ix 



xsm 



o 



knx 



kir 

2 



-xcos 



2 
knx 



kirx 



dx 



7 cos(kir) + -s— ^ 

K7T K Z 7T Z 



dx + l 811) 

1 r 1 2 /jfeTra;', 

+ / — - cos ^— I dx 

/C7T V 2 

fc7T N 







2 ( ki\x 

— cos 

kir V 2 



sin 



4 sin 



2fc7r( — 1)' 



2^2 



fc 2 7T 
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Conclui-se que o desenvolvimento de f em série de senos é dado pela expressão 

e . . ^4sin(^)-2fc7r(-l) fc . / kirx 

/(*) = £ — ^ **{-t- 

fe=i v 

(b) O desenvolvimento de f em série de senos no intervalo [0, 2n] é da forma 

f( x ) = X^fcsinf — \ . 

Os coeficientes bk são dados pela fórmula 

2 f 27T ., N . fkx 



bu = — / f(x) sin — dx 
2tt 7o V 2 ) 

i r 2 - . (kx\ 

= — cos i sm — dx. 

* Jo V 2 y 

Para calcular a primitiva anterior pode integrar-se duas vezes por partes e obtém-se 

k 2 \ f ( kx\ f kx\ k fkx 

1 — I / cos x sm I — I dx = sm x sm I — I + — cos x cos 

Assim, para k ^ 2, tem-se 

1 4 / . . ( kx\ k ( kx \ \ 



bu = — ■ ft sin x sin — M — cos x cos , 

tt 4-k 2 \ V 2 ) 2 \2 

2((-l) fc -l)fc 



^(4-fc 2 ) 
e, para k = 2, tem-se 



1 P w 1 /" 27r , . 

02 = — / cosa; sin x dx = — / sin(2x) dx = 0. 

71" Jo 2tt 7o 

Conclui-se que o desenvolvimento de f em série de senos é dado pela expressão 



«■>- E 2 ^Mf) 



fc=l,fc^2 

D 



(3) Determine o desenvolvimento em série de cosenos das seguintes funções: 

(a) f(x) = x 2 para x e [0, 1]; 

(b) /(a;) = e 2 * para x e [0,2tt]. 

Resolução: 

(a) O desenvolvimento de f em série de cosenos no intervalo [0,1] é da forma 

oo 

f( x ) = y + £ ak cos ( fcvra; )- 

k=l 
Os coeficientes ak são dados pela fórmula 

2 f 1 
ãk = - I f(x) cos(k7:x) dx. 
1 Jo 

Donde 

f 1 i 2 

ao = 2 / x dx = -, 
Jo 3 
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e, para k > 1, 



a/j = 2 x cos(knx) dx 

/o 

., isllll irn.r i 

= 2\ x 



kn 



o Jo 



1 sin(/í7rx) 



Ã'7T 



xsm(kirx) dx 
i 



4 í x cos(A;7rx 

fc7T l 

4(-l) fe 



kn 



1 COs(/C7Tx) 

+ / ^ dx 



n 2 k 2 ' 
Conclui-se que o desenvolvimento de f em série de cosenos no intervalo [0,1] é dado 
por 



/w4 + Ê íHI ' 



3 ' ^ 7T 2 fc 2 

k = \ 



cos(knx). 



(b) O desenvolvimento de f em série de cosenos no intervalo [0, 2ir] é da forma 

oo 
+ 22 &k COS 



m 



ao ^- / kx 

~2 + 

k=\ 



Os coeficientes a& são dados pela fórmula 

at = — / f(x) cos ( — ) dx. 
2tt 7o V 2 



Donde 



e, para k > 1, 



a = 



7T J 



2- 



e 2x dx = 



«A- 



1 

7T J0 

= -Re 

7T 



2- 



e 2:E cos 



e ivr -1 
fcx\ 



2- 



e 2 * +i " dx 



o 



i / e( 2 +f*> 
= —R e i. — 

7T \ 2+ |i 



1 / p A7T-\-k7TÍ 

= -Re r - 

7T \ 2+£í 



2ttN 



8((-l) 



k„4ir 



1 



vr(16 + fc 2 ) 
Conclui-se que o desenvolvimento de f em série de cosenos é dado pela expressão 



f(x) 



e An_ 1 ~ g ((_!) 



fe p 47T 



47T 



fc=l 



1) /jfez\ 

7TT COS . 

vr(16 + fc 2 ) \2 ) 



D 



AMIV - FICHA RESOLVIDA 6 5 

(4) Recorrendo ao método de separação de variáveis determine uma solução do seguinte 
problema de valor na fronteira para a equação do calor: 

du d 2 u 

dt ~ dx 2 

u(0,x) = sin(3a;) — | sin(8x) 
u(t,0) =u(í,7r) = 

para 0<x<7reí>0. 

Resolução: Começa-se por procurar soluções da equação diferencial parcial da forma 
u(t,x) = T(t)X(x). Substituindo na equação tem-se 

§- (T(t)X(x)) = |^ (T(t)X(x)) 
ot ox z 

<=> T'(t)X(x) = T(t)X"(x) 
T'{t) X"{x) 



T(t) " X(x) 



paraT(t),X(x) ^ 



riu \ i ví \ P ara a lg um k £ 



T'(t) , X"(x) 

Wf = k= im para algum ke 

T'(t) = kT(t) 
X"(x) = kX(x) 

Tendo em conta as condições na fronteira 

u(t,0) = u(t,ir) = ^^ T(t)X(0) = T(t)X(n) = 

T(t) = Ví ou X(0) = X(tt) = 0, 

conclui-se que para obter soluções não identicamente nulas tem de ser X (0) = X(tt) = 0. 
Por sua vez isto implica que a constante k acima tem de ser negativa (se k > a única 
solução da equação X"(x) — kX(x) = que verifica X(0) = X(ir) = é a solução 
identicamente nula). 
Portanto 

X(x) = ci cos(V-kx) + C2 sm(V-kx), 
e substituindo nas condições fronteira obtém-se 

X(0) = _ íci = í ci = 

X(tt) = | C2 sin(v / — kir) = \ C2 = ou sin(v / — kn) = 0. 

Obtêm-se assim soluções não nulas para a equação diferencial para X e para a condição 
fronteira quando 

I 9 

V —kir = n-K <í=> k = —n para n = 1,2, .. . 

e para cada um destes valores de k obtêm-se como soluções do sistema acima múltiplos 
reais de 

X(x) = sin(ra) e T(t) = e~ nH . 

Isto é, para cada n = 1,2,... obtém-se a seguinte solução da equação diferencial parcial 
satisfazendo a condição na fronteira: 

u n (t,x) = sin(na;)e _n *. 

Finalmente, determina-se coeficientes d„eK tais que a condição inicial seja satisfeita por 



u(t,x) = y d n u n (t,x). 



n=\ 



Tem-se 
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u(0,x) = sin(3x) — -sin(8x) 



oo 1 

y d n u n (0, x) = sin(3x) — - sin(8a;) 

re=l 

oo 1 

y d n sin(nx) • 1 = sin(3a;) sin(8x) 



re=l 



portanto efe = 1, efe = — \ e d n = para n ^ 3,8. Conclui-se que a solução do problema 
do enunciado é 



u(t,x) = sin(3x)e" 9í - - sin(8a;)e" 64í . 



D 



(5) Seja c um parâmetro real positivo. Recorrendo ao método de separação de variáveis 
determine uma solução do seguinte problema para a equação das ondas: 

d 2 u _ 2d 2 u 

at 2 c dx 2 
it (0, x) = cosx 

§1(0,*) = o 

u(t,Q) = u(t,2ir) = 1 

para 0<x<27reí>0 (verificando a equação diferencial para < x < 2ir). 

Sugestão: Comece por determinar uma solução estacionária (isto 
é da forma u(t, x) = v(x)) da equação diferencial que satisfaça as 
condições na fronteira para x = e x = 2ir. 
Pode também aproveitar o resultado da alínea 2(b). 

Resolução: Começa-se por determinar uma solução estacionária da equação diferencial 
parcial que satisfaz a condição na fronteira. Substituindo u(t, x) = v(x) em 

d 2 u 2 d 2 u 

dt 2 ~ dx 2 

tt(í,0) = u(t,2ir) = 1 

obtém-se 

g?(v(x)) = c 2 £ T (v(x)) f v(x) = ax + b .,m = i 

v(0) = v(2n) = 1 1 u(0) = v(l) = 1 ^^ V{X> i ' 

Voltando ao problema inicial e escrevendo 

u(t, X) = v(x) + Uh(t, X) = 1 + Uh(t, x), 

tem-se que Uh(t,x) é uma solução do seguinte problema com condições na fronteira ho- 
mogéneas: 

d 2 u h = c 2 d 2 u h 
dt 2 dx 2 

Uh(0,x) = cosx — 1 

u h (í,0) = u fc (í,27r) = 
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Para resolver este problema, começa-se por procurar soluções da equação diferencial parcial 
e das condições fronteira da forma Uh(t,x) = T(t)X(x). Substituindo na equação tem-se 

^(T(t)X(x)) = c^ 2 (T(t)X(x)) 
<=► T"{t)X{x) = c 2 T(t)X"(x) 

~ m=^» araTit) ' x(x) *° 

T\t) 2 X "(x) 

<^> -f-( = k = c 2 —^- para algum keR 
T(t) X(x) 

T'(t) = kT(t) . , _ 

X"(x) = $X{x) ^ dlgUm k G R - 

Tendo em conta as condições na fronteira 

u fc (í,0) = u h (t,2ir) = ^> T(t)X(0) = T(t)X(2ir) = 

T(t) = Vi ou X{0) = X(2tt) = 0. 

conclui-se que para obter soluções não identicamente nulas tem de ser X(0) = X(2tt) = 0. 

Por sua vez isto implica que a constante k acima tem de ser negativa (se k > a única 

solução da equação X"(x) — ^(X(x) = que verifica X(0) = X(2tt) = é a solução 

identicamente nula). 

Portanto 

fV^k \ . íV^k 

X (x) = c\ cos I x I + C2 sin 

e substituindo nas condições fronteira obtém-se 

X(0) = í ci = _ 



X(2tt) = I c 2 sin V^-2tt = 0. 




Obtêm-se assim soluções não nulas para a equação diferencial para X e para a condição 
fronteira quando 

2\/—kTT = ncTT <í=> k = — para n= 1,2,... 

e para cada um destes valores de k obtêm-se como soluções da equação para X(x) 

X(x) = c 2 sin (^—j . 
Quanto a equação para T(t) 

Tl C ( Tld \ ( Tict \ 

T"(t) = -—T(t) <=> T(t) = c 3 cos (^— J + c 4 sin (^— J 
No entanto, a condição inicial 

d ^(0,x) = T'(0)X(x) = 

implica, para soluções não identicamente nulas, que seja T"(0) = 0, ou seja, C4 = 0. Assim 

' nct s 



T(t) = c 3 cos , 9 
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Isto é, para cada n = 1,2,... obtém-se a seguinte solução da equação diferencial parcial 
satisfazendo a condição na fronteira e a condição inicial respeitante a derivada em ordem 
a t: 



. . /nx\ 

u n (t,x) = sm I — I cos 



nct\ 



2 ) 

Finalmente, determina-se coeficientes 4éR tais que a condição inicial não homogénea 
seja satisfeita por 

OO 

u h (t,x) = y^d n u n (t,x). 

71=1 

Tem-se 



Uh(0,x) = cosa; — 1 

oo 

y d n u n (0, x) = cos x — 1 



n=\ 

oo 



71=1 



V 2 



y d n sin í — ) • 1 = cos x — 1 . 



Portanto os d n 's são os coeficientes do desenvolvimento da função cos x — 1 em série de 
senos no intervalo [0,2tt]. Na alínea 2(b) foi já calculado o desenvolvimento de cosa; em 
série de senos neste intervalo pelo que resta fazer o mesmo para a função constante igual 
a -1: 



2~ 



ir(-i)»n*=--n 1, - !(l41 '"" 

2tt 7o ^ 2 / nvr V 2 / 



n-K 



donde se conclui 

oo 



2((-l) ra -l) fnx\ 

sm — 

nvr V 2 I 

Tt=l 



e portanto, tendo em conta que para n par os coeficientes das séries de senos das funções 
cosa; e —1 se anulam, tem-se 

Ef An 4 \ fnx 
— —r. ?rr H sm — 
\7r(4-n 2 ) mrj \ 2 

n=í,n ímpar 

isto é, 

An 4 



se n impar 
Conclui-se que 



d n = i T(4- ra 2) nvr 

sen par. 



^— v / 4n 4 \ /na;\ f nct\ 
Uhit, x) = > — — — =- H sm — cos 

n=l, n impar 



impar 

e que a solução do problema do enunciado é 



. . v-^ ( 4n 4 \ . /ra\ / nci 
u\t,x) = 1+ > — — ; =- H sm — cos — 

n=í,n impar 



D 
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(6) Recorrendo ao método de separação de variáveis determine uma solução para o 
seguinte problema de valor na fronteira: 

d 2 u i <9 2 ? 



dx 2 ' dy 2 

u(x, 0) = u(x, 1) = u(0, y) = 

u(í,y) = y 



'ii ( \ — i ví \ P ara a lg um k ^ 



para 0<x<le0<y<l (verificando a equação diferencial para < x < 1 e 
Q<y<l). 

Resolução: Começa-se por procurar soluções da equação diferencial parcial da forma 
u(x,y) = X(x)Y(y). Substituindo na equação tem-se 

J^ (X(x)Y(y)) + ^ (X(x)Y(y)) = X{x)Y{y) 

^ X"(x)Y(y) + X(x)Y"(y) = X{x)Y{y) 
X"(x) Y"(v) 

^^ —ff -l = k = -^tt^v P ara al êum k£R 
X(x) Y(y) 

X"(x) = (k + l)X(x) 
Y"(y) = -kY(y) 

Tendo em conta as condições na fronteira 

u(x, 0) = u(x, 1) = ^^ X(x)Y(0) = X(x)Y(l) = 

x( x ) = ovx ou y(o) = r(i) = o 

conclui-se que para obter soluções não identicamente nulas tem de ser Y(0) = Y(l) = 0. 
Por sua vez isto implica que a constante k acima tem de ser positiva (se k < a única 
solução da equação Y"(x) + kY(y) = que verifica Y(0) = Y(l) = é a solução 
identicamente nula). 
Portanto 

Y(y) = ci cos(Vky) + C2 sin(v ky) 
e substituindo nas condições fronteira obtém-se 

y(o) = o _ f ci = o í ci = o 

Y(l) = ^^ \ c 2 sm(Vk) = ^^ \ c 2 = ou sm(Vk) = 

Obtêm-se assim soluções não nulas para a equação diferencial para X e para a condição 
fronteira quando 

Vk = rnr <í=> k = n 2 ir 2 para n = 1,2,... 

e para cada um destes valores de k obtêm-se como soluções da equação para Y{y) 

Y{y) = C2 sm(niry) 

Quanto a equação para X(x) 

X"(x) = (l + nV)I(s) 

<í=^» X(x) = C3 cosh( v 1 + n 2 ir 2 x) + C4 sinh( y 1 + n 2r K 2 x) 

No entanto, a condição 

u(0,y) = X(0)Y(y) = 
implica, para soluções não identicamente nulas, que seja X(0) = 0, ou seja, C3 = 0. Assim 



X(x) = C4sinh(v 1 + n 2 TT 2 x). 
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Isto é, para cada n = 1,2,... obtém-se a seguinte solução da equação diferencial parcial 
satisfazendo as condições na fronteira homogéneas: 

u n (x, y) = sinh( v 1 + n 2 n 2 x) sin(mry) 

Finalmente, determina-se coeficientes d n 6 R tais que a condição na fronteira não ho- 
mogénea seja satisfeita por 

oo 

u(x,y) = y]dnU n (x,y). 

ra=l 

Tem-se 

u(l,y) = y 

oo 

^^ ^d ra u„(l,i/) = y 

Tl=l 

OO 

<í=> Y^ d n sinh( v 1 + n 2 7r 2 ) sin(n7ry) = y. 

n=\ 

Portanto d n sinh(v / l + n 2 ir 2 ) são os coeficientes do desenvolvimento de y em série de 
senos no intervalo [0,1]. Donde, 

/ 2 f 1 2(— l) n+1 

d ra sinh(v 1 + n 2 7r 2 ) = - / ysmímry) dy = 

1 Jo rnr 

e portanto 

2(_i)«+i 



n7r sinh(v / l + n 2 7r 2 ) 
Conclui-se que a solução do problema do enunciado é 

°^ 2(— l) n+1 / 

u(x, y) = 2, — sinh(v 1 + n 2 Tr 2 x) sm(mry). 

~ l mrsinh.(v 1 + n 2 7r 2 ) 



D 



